Intersubband Electron Interaction in 1D-2D Junctions 
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We have shown that the electron transport through junctions of one-dimensional and two- 
dimensional systems, as well as through quantum point contacts, is considerably affected by the 
interaction of electrons of different subbands. The interaction mechanism is caused by Friedel oscil- 
lations, which are produced by electrons of the closed subbands even in smooth junctions. Because 
of the interaction with these oscillations, electrons of the open subbands experience a backscatter- 
ing. The electron reflection coefficient, which describes the backscattering, has a sharp peak at the 
energy equal to the Fermi energy and may be as high as about 0.1. This result allows one to explain 
a number of available experimental facts. 



Quantum wires and quantum point contacts are of 
prime interest as model systems for studying the effects 
of electron-electron interaction, which is known to play 
an important role in one-dimensional systems, giving rise 
to a correlated state. Strong evidence in favor of the 
Luttinger liquid was obtained from the studies of tunnel- 
coupled quantum wires P, 0|. However, experiments 
with isolated quantum wires do not directly testify to the 
presence of a Luttinger liquid but reveal many transport 
features that have not yet found any adequate explana- 
tion. Presumably, these features are related not only to 
electron-electron interaction in a wire, but also to the 
fact that the wire is connected with electron reservoirs. 
The most prominent and most discussed feature is the 
anomalous conductance plateau at 0.7 x 2e /h$. Other 
experimental facts can be classed into three groups: 

(i) Experiments testifying to the electron localization 
above the potential barrier formed by a smooth (on the 
Fermi wavelength scale) electrode potential [1, @| ■ Such 
a localization is supposed to interpret the 0.7 x 2e 2 /h 
anomaly in terms of the Kondo effect [1, [f|. The lo- 
calization mechanism remains unknown. Numerical cal- 
culations Q performed to justify the spin localization 
assumption have little force, because they use geomet- 
ric dimensions of the contact that are close to the Fermi 
wavelength or even smaller. 

(ii) Studies of the nonlinear conductance at a small 
(compared to the intersubband and Fermi energies) volt- 
age applied along the quantum wire. An increase in the 
height of the differential conductance plateau is observed, 
whereas the ballistic conductance theory [H, Q predicts 
its decrease with voltage and the inclusion of the electron- 
electron interaction via the selfconsistent field does not 
qualitatively change this conclusion Moreover, in 
the experiment, the first conductance plateau rises to a 
level even higher than 2e 2 /h [HI, 12|, which points to 
the appearance of an additional transport channel, e.g., 
through higher size-quantization subbands. 

(iii) Observation of a specificc scattering in the regions 
of sufficiently smooth transitions between one- and two- 
dimensional electron systems (1D-2D junctions). The 
scattering manifests itself as an effective resistance (esti- 



mated as ~ 0.1/t/2e 2 ) connected in series with the quan- 
tum point contact [12| or as the effect of the potential 
profile of the junction on the structure of the 0.7 • 2e 2 jh 
anomaly [HI]. 

The present paper shows that these experiments can 
be explained (at least qualitatively) if one takes into ac- 
count the interaction between electrons of different sub- 
bands in the junction between the ID and 2D parts of 
the structure and, primarily, the interaction of the elec- 
trons passing through the junction with the electrons of 
the closed subbands. The interaction mechanism is re- 
lated to the Friedel oscillations of electron density, which 
occur in the junction because of the reflection of higher 
subband electrons not passing through the contact. The 
physical picture is as follows: numerous electrons in the 
reservoirs collide with the contact, but only a small num- 
ber of them (electrons belonging to the open subbands) 
can pass through the contact. All the other electrons 
are backscattered, causing Friedel oscillations of electron 
density. These oscillations evidently have different phases 
in different subbands. The phases depend on the form 
of the junction, but, since this form is described by a 
regular function, the summation over the subbands does 
not lead to the disappearance of oscillations. Away from 
the contact, the oscillations have a clearly pronounced 
component with a wave vector 2kF (kp is the Fermi 
wave vector in the 2D reservoir). This assumption is 
supported by the experiment, which reveals the oscilla- 
tory structure of the electron density distribution at a 
large distance from the contact with the use of a probe 
microscopy technique FJ]. Our calculations show that 



the interaction of the electrons passing through the con- 
tact with the Friedel oscillations leads to a fairly strong 
backscattering. With this fact taken into account, it is 
possible basically to explain the experiments mentioned 
above. 

Let us consider a 1D-2D junction in the form of a strip 
(Fig. 1) whose width d(x) monotonically increases along 
the x axis from d > kp 1 at x — to D d at x — > oo. 
The characteristic expansion radius R considerably ex- 



ceeds both d and k. 



The problem consists in the 



evaluation of the reflection coefficient of electrons that 
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FIG. 1: Electron density oscillations in a smooth 1D-2D junc- 
tion: (a) the strip width variation in the junction, (b) the size 
quantization subbands, and (c) the electron density oscilla- 
tions in the closed subbands (n = 2, 3, 4 . . . ). 



are incident on the contact in the open subband and re- 
flected as a result of their interaction with the Friedel os- 
cillations caused by the electrons of the closed subbands. 
For calculations, we use the Born approximation, which 
is justified if the reflection coefficient is small. In the 
zero-order approximation with respect to the interaction, 
the wavefunctions can be determined in the framework 
of the standard adiabatic approximation (see, e.g., [Hj]). 
For the closed subbands, we have 



To determine the potential of the perturbation causing 
the transition from the |n, k) state to the |m, k') state, we 
use the HartreeFock approximation. Taking into account 
that the 1D-2D junction is smooth, we reduce the prob- 
lem to effective one-dimensional equations by integrating 
the HartreeFock equations with respect to the transverse 
coordinates. As a result, we obtain the following expres- 
sion for the reflection coefficient for the electrons in the 
open subband (i.e., for the (1, k) — > (1, —k) transition): 



where %pk is the ^-dependent part of the function 
^i,k{x, y) in Eq.([T]). The perturbation potential contains 
the Hartree and exchange components: V = V H + V exc . 
The Hartree potential is the function 



V H (X)- 



£ 

n=2 



dx'Vf n (x,x')p n (x')~ dx'V" (x,x')p (x') , 



where n(x) is the electron density in the nth subband and 
po is the positive background charge density. The per- 
turbation caused by the exchange interaction is described 
by the operator 



V exc ^{x) = - 



E 

71 = 2 



dx'Vr n C (x^')Pn(x,x')^(x'), 



where Nf is the index of the upper subband filled in the 
2D part of the system and p n (x,x') is the density ma- 
trix, or, more precisely, its perturbation caused by the 
1D-2D junction. The effective potentials V^ m (x, x') and 
V^^(x,x') for the direct and exchange interactions be- 
tween the electrons belonging to the nth and mth sub- 
bands have the form 

V n H m (x, x') = [ d(X) dy dyVfo r') X 2 n (y) X 2 m (y') , 
Jo Jo 



V H (x,x') = d(x')- 1 [ diX) dy [ d(X dy'V(r,r')xl(y) , 
Jo Jo 



k n (x) 



dx'k n (x') - — 



where n = 2, 3, . . . is the number of a subband; k n {x) 
is the wave vector of the longitudinal motion; k = 
lima— ,oo k n (x); Xn(y) is a transverse wavefunction; and 
a n {k) is the turning point. For simplicity, we assume 
that only one subband is open. For this subband the 
wave function is 



*i,fe(x,y) ~ Xi{y)\ 



k 



kx(x) 



exp 



i I dx'k n (x') 
o 



(1) 



d(x) rd(x') 

dy \ dy'V{r, r') x 
o Jo 

(y)xm{y)xn(y')xm{y') ■ 



Here, V(r, r') is the pair interaction potential, which 
is determined by the Coulomb interaction screened by 
metal electrodes, if any, and by two-dimensional elec- 
trons. 

The calculation of the reflection coefficient can be sim- 
plified by taking into account the actual structure of the 
spatial distribution of electron density. Two density com- 
ponents are present: p(x) ~ p(x) + p(x), where the first 
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component smoothly varies on the kp 1 scale and the sec- 
ond component oscillates with a wave vector of ~ 2kp, 
with the oscillation amplitude smoothly (approximately 
as a; -3 / 2 ) decaying toward the depth of the 2D region. 
It is of special interest to consider the far zone (x > R) 
where the wave vector of oscillations is close to 2kp, be- 
cause these oscillations most efficiently scatter the elec- 
trons with Fermi energy in the backward direction. The 
contribution of the smooth component p(x) is small in 
terms of the parameter (Rkp)^ 1 -C 1. The contribu- 
tion of the oscillating component p{x) in the near zone is 
unimportant, because, here, the oscillation period is no- 
ticeably greater than the electron wavelength in the open 
subband. Thus, the problem can be simplified for the re- 
gion x > R with allowance for the fact that the density 
oscillations of interest with the wave vector ~ 2kp are 
produced by the electrons belonging to the lower sub- 
bands (n <C Dkp/n) and characterized by a longitudinal 
momentum close to kp . 

For such electrons, the interaction potentials 
Vj^ m (x,x') and V 7 ^(x,x') are simplified if the ef- 
fective pair interaction radius a is small compared to 
the wavelength in the transverse direction. In reality, 
this condition is satisfied in the far zone, because, in 
this zone, a is on the order of the Bohr radius as and 
the transverse wavelength is on the order of D/n. In 
this case, the difference between the potentials V^ m and 
V^ X m i as well as their dependence on the band indices, 
is insignificant. In addition, at distances greater than a, 
the interaction potential can be assumed to depend on 



the coordinate difference 



Thus, we obtain 



V H « V exc « V(x - x 1 ) = —U ( ( 2 ) 

e D \ a J 

where eo is the permittivity of the semiconductor and 
U(x) is the dimensionless potential, which depends on 
the screening in the system. In the case of the Coulomb 
interaction, the effect of screening by 2D electrons can 
be taken into account [1 61 ] by using the permittivity of 
the 2D electron gas e(q) in the random phase approxi- 



mation 17 1 



In this case. 
U(x) = i 



dq 



\Q\<Q) 



The oscillating electron density component p{x) is de- 
termined by the summation over all of the closed sub- 
bands. If the electron-electron interaction is ignored, the 
oscillation amplitude caused by an individual subband 
decreases with distance as 1/x. The inclusion of interac- 
tion leads to a slower decrease [H, [l9[ of the amplitude. 
Because of the variance in the oscillation phases of dif- 
ferent subbands, the total amplitude proves to decrease 
with distance faster than 1/x. For example, in the case 
of a at edge of the 2D gas, the amplitude asymptotically 
decreases as x _3 / 2 [l(|. The presence of the 1D-2D junc- 
tion leads to a greater phase variance and, hence, to a 
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FIG. 2: Phase vs. the band index (the solid curve) for the 
1D-2D junction and (the dashed line) for the straight edge of 
the 2D electron gas. 



decrease in the Friedel oscillation amplitude. In the re- 
gion x > R, for the electrons of the closed subbands, the 
phase at the Fermi level can be represented as 

$(7, x) ~ k F x\J\ -7 2 - $(7, k F ) , 

where 7 = nn/ (kpD). If kpD 3> 1, one can assume that 
7 is a continuous quantity and replace the summation 
over n by integration with respect to 7. An important 
role is played by the second term $(7, kp) associated 
with the presence of the 1D-2D junction. In the absence 
of the junction, i.e., for a straight edge of the 2D region, 
$(7, kp) — 0. The specific form of the function $(7, kp) 
is determined by the form of the junction, but the general 
property of 1D-2D junctions is a sharp increase of the 
phase at 7 — > 0, as it is shown in Fig. 2. 

To make the consideration more specific, let us con- 
sider the case specified by 



d{x) 



d/VT 



(x/R) 2 , 



x < rVT 
x > rVT 



6 2 



S 2 , 



(3) 



where S = d/D <C 1. The characteristic feature of this 
d(x) dependence is the presence of almost straight edges 
away from the transition region. Let us denote £ = x/R. 
In the region £ > \/ 1 — S 2 , 

S(7,aO~Jtt* [£-l + <5 2 /(37 2 )] 

Calculating the density with the use of asymptotic ex- 
pansions in Xp = 2Rkp 3> 1, we arrive at the following 
result: 



p(x) 



2k F D riTcxp-\ F 6y/2(£-l)/3 



7T AttR V 2\ F i - 1 

H( "e X p{z[A F (£-l)-7r/4]} 



One can see that, unlike the case of the straight edge, 
in the 1D-2D junction, the density oscillations decay 
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exponentially with a characteristic decay length of ~ 
(D /kpd) 2 1 R. For D ^> R, this length is much greater 
than R and, hence, kp 1 . Therefore, such a decay of os- 
cillations is insignificant. 

The reflection coefficient, correct to the phase factor, 
is expressed as 



akf 



it\/tt ask 



U2ka~ 



Uo 
2 



F 



k 

kp 



A[*,[/] , (4) 



where the Fourier transform U q of the potential U(t), 
which depends on the dimensionless coordinate t = x/a 
is introduced; Uika and U§ are U q at q = 2ka and q = 0. 
Note that, although the potentials V H and V exc (see Eq. 
(J2j) ) depend on D, the final expression for ru is free of this 
dependence because of the summation over the subbands. 
The function F(k/kp) has the form 



F(z) 



r cos(2y + c-ir/4) 2izv „ 2SJXf(v+c/2)/3 
J (y + c/2)3/ 2 



where c ~ 1 is a constant arising because of the limita- 
tion imposed by the applicability of the asymptotic for- 
mulas for Friedel oscillations. In Eq. ((4]), the term U2ka is 
caused by the Hartree interaction, and the term Uq is as- 
sociated with the exchange interaction. The combination 
\2Uika — Uq] appears in a standard way in the theory of 
electron scattering by Friedel oscillations 1(| ID, 2^, 21 1 . 
The last term in Eq. g]), 

A[ZiU] = r dy C ™^ +C -Z^ e^W)/s x 
Jo {y + c/2y/ 2 



dtU(t), 

2y/ak F 

is also associated with the exchange interaction. It arises 
because of the difference in the direct and exchange in- 
tersubband interactions in the 1D-2D junction. Within 
the part of the junction, where the electron density of 
closed subbands is absent, the exchange interaction be- 
tween the closed subband electrons and the electrons of 
the open subband is also absent, while the direct inter- 
action extends over a distance of about a. For o/cf < 1, 
the term A[z, U] vanishes. In any case, A[z, U] does not 
noticeably affect the reflection coefficient, so that \rk \ can 
be estimated by the first term of Eq. ([4]). 

The dependence of the reflection coefficient on the 
wave vector k is mainly determined by the function 
(fcp /k)F(k/kp) plotted in Fig. 3. The characteristic 
feature of this function is the sharp peak at k = kp, 
near which the function follows a root dependence on 
\k — kp\. The factor \U2ka — Uo/2] varies more smoothly. 
If U(t) = U a exp(-\t\), we have U = 2U a and U 2 ka = 
2U a /(l + 4k 2 a 2 ). In the case of a screened Coulomb in- 
teraction, Uq ~ it/2 and U2ka — 7r /[2(l + fcas)]- Hence, 
the dependence of | | on A; is approximately identical to 
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FIG. 3: Function (ItF /k)F(k/kp) determining the depen- 
dence of the reflection coefficient on the electron energy; 
\f — 30, 5 — 10~ 3 , c = 0.5. 



the dependence shown in Fig. 3. The value of \rk\ for 
actual values of the parameters kp and a can be on the 
order of several tenths, and |rfc| 2 ~ 0.1. The contribution 
of the exchange interaction is predominant. 

The results obtained above qualitatively hold for 1D- 
2D junctions of other forms. For example, in the case of a 
smoother junction with d(x) = D — (D~d)exp[—(x/R) 2 ], 
the phase $(7, kp) has a root singularity for small val- 
ues of 7: &(-f,kp) ~ RkF^ — ci-y/j) yl — 7 2 , where 
ci ~ 5/4. In this case, a length parameter appears, 
l c — 4(cii?) 4 /cp, which is considerably greater than R. 
Since the phase variance is stronger than that in the case 
considered above, i.e., in the case of d(x) given by Eq. ([3]), 
the Friedel oscillation amplitude decreases but its depen- 
dence on distance weakens within the length l c : p ~ x . 
For x l c , the dependence p ~ x~ 3 l 2 holds. As a re- 
sult, the electron reflection coefficient characterizing the 
scattering in the far zone decreases approximately by a 
factor of ~ it 2 (ci\f)~ 2 \n(aAp). However, at the same 
time, the size of the near zone increases and the role of 
the scattering processes in it becomes more significant. 
This case requires special consideration. 

Thus, the interaction between electrons of different 
subbands in 1D-2D junctions can be sufficiently strong 
to affect the electron transport in quantum point con- 
tacts and quantum wires. The effect of this interaction 
is as follows: 

(i) The backscattering of electrons of open subbands 
leads to a decrease in conductance, which explains the 
presence of a specific resistance observed experimentally 
for 1D-2D junctions 



12] 



The estimate obtained above, 
\fk\ 2 ~ 0.1, agrees well with these experiments; 

(ii) In a quantum wire connecting two electron reser- 
voirs, the electron scattering occurs in two opposite junc- 
tions, which results in the appearance of quasibound 
states and enhances the effect of the scattering mech- 
anism under discussion on the conductance; 

(iii) Since the scattering by Friedel oscillations de- 
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creases with an increase in the electron energy for k > hp, 
the differential conductance of the contact may increase 
with increasing bias voltage. In principle, the differen- 
tial conductance may even exceed the value of 2e 2 //i, 
because the intersubband interaction creates a possibil- 
ity for the transitions from the closed subbands to the 
open ones (even for smooth fD-2D junctions). Indeed, 
the matrix element of the electron transition from the 
nth subband to the first one due to the interaction with 
the electrons of the mth subband is determined by the 
terms of the form (xi(y)Xm(y)Xn(y)) 7^ 0. If, for exam- 
ple, Xn(y) ~ sin(7rm/ / ' D) , such transitions are possible 
for all of the odd values of n. 
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